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Abstract. In this paper, using a standard quantum chemistry approach based on 
pseudopotentials for atomic core representation, Gaussian basis sets, and effective 
core polarization potentials, we investigate the electronic properties of the MgH+ ion. 
We first determine potential energy curves for several states using different basis sets 
and discuss their predicted accuracy by comparing our values of the well depths and 
position with other available results. We then calculate permanent and transition 
dipole moments for several transitions. Finally for the first time, we calculate the 
static dipole polarizability of MgH+ as function of the interatomic distance. This 
study represents the first step towards the modeling of collisions between trapped cold 
Mg+ ions and H2 molecules. 
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1. Introduction 

The reactivity of molecular systems in the gas phase at low temperature (as low as 
about 10 K) is of crucial importance for the evolution of cold natural environments 
in the interstellar medium. In the laboratory, amazing progresses have been achieved 
in the obtention of gaseous molecular samples at even lower temperatures, from the 
cold regime (a few K or less) down to the ultracold regime (a few mK and below). 
It is now possible to study chemical reactivity in the Kelvin range down to 6 orders of 
magnitude in temperature, with a wide variety of molecular systems [1]. In this context, 
molecular ions are appealing as they can be easily manipulated by electromagnetic fields 
to guide and confine them in traps, until the formation of ion crystals [2]. Nowadays, 
cold molecular ions can be produced by sympathetic cooling with laser-cooled ions [3] 
or by collisions with a buffer gas of helium of hydrogen [1]. Interactions between cold 
atoms and ions are also relevant for open question related to molecule formation in 
Bose- Einstein condensates [5] and to quantum information 

The first evidence for the formation of molecular ions in a laser-cooled ion trap 
has been reported by Baba and Waki [7], by introducing air in a Mg^ trap. A more 
controlled experiment has been achieved by Drewsen and coworkers [8] , who introduced 
thermal H2 or D2 gas in a laser-cooled Mg"*" trap yielding trapped MgH"*" or MgD"*" 
ions. Other diatomic or triatomic molecular species have been produced, which are for 
instance reviewed in ref. [9] . The main challenge is the detection of the formed molecular 
species, which is most often indirect, but non-destructive. Indeed, the molecular ions 
are not sensitive to the cooling laser, so they do not contribute to the fluorescence of 
the trap. By controlling the gaseous species which is introduced in the cold ion trap, 
and knowing the open reactive channels, one can observe the absence of fluorescence 
at specific locations of the ion crystals, assigned to the presence of molecular ions. 
The changes in the spatial structure and the excitation by an oscillating voltage [10] 
of the kinetic modes of the dual-species ion crystal can also be detected, and modeled 
considering the masses of the expected molecular ions. A more direct detection method 
- probably destructive- is needed, relying on the spectroscopic characteristics of the 
molecular ions. This would also allow to address the other major challenge of such 
experiments, i.e. the knowledge and the control of the internal energy of the formed 
molecular ions after the reactive process. 

In this perspective, the MgH"*" ion represents a prototype system, which is 
extensively studied by the Aarhus group. Bertelsen et al [llj observed for the first 
time the photodissociation of trapped MgH^ ions embedded in a Mg"*" ion crystal, using 
a two-photon resonant scheme. The simulation of the MgH"*" photodissociation and the 
possibility to control the branching ratio has been investigated by J0rgensen et al |12j . 
The isotope effects in the reaction of excited Mg"*" ions with either H2 of HD molecules 
have been observed [13] through a factor of 5 difference in the branching ratios for MgH+ 
and MgD+ ions. 

It is also worthwhile to mention that prior to the above works, the effective two- 
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electron nature of the MgH"*" molecule made it appropriate as a prototype for accurate 
quantum chemistry calculations at various stages of development of computing facilities 
[m [iSl [ini [IT]. The MgH"*" system also attracted attention in various situations of 
collision physics. Charge transfer in collision between H"*" and Mg have been investigated 
in the keV energy range [iHl [17], and related cross sections have been calculated |T5] . 
Broadening of the Mg"*" resonance line by collision with atomic hydrogen has been 
investigated by Monteiro et al [19j. Dalleska et al investigated the energy dependence 
of the Mg++H2 reactive collision and estimated the MgH"*" dissociation energy |20j . 
In most of these works, spectroscopic constant for various states of MgH"*" have been 
reported. Nevertheless, to our knowledge, there is only one high-resolution spectroscopic 
study of MgH+, performed by Balfour in 1972 [21j. The rotational analysis of the 
A^E+ X^S+ and B^U transitions of several isopotomers (2''MgH+, ^^MgH+, 

^^MgH+) enabled the author to derive the low part of the potential curves for the X and 
the A states, and some spectroscopic constants for the B state. According to ref . [T^ . 
another spectroscopic analysis has been performed by Numrich [22] , but is not available 
to us. 

In this work, we investigate electronic properties of MgH+ ion as a first step towards 
the modeling of collisions between trapped cold Mg"*" ions and H2 molecules. Following 
our work on effective two-electron diatomic molecules like alkali dimers [23l [2H [25] 
and alkali hydrides [26], we determine potential curves for various symmetries up to 
the Mg"'"(5s)-fH(n=l) asymptote, static dipole polarizabilities, and permanent and 
transition dipole moments for the main transitions involving the lowest electronic states 
relevant for photodissociation. In several occurrences, such results are given for the first 
time. We will often use atomic units for distances (ao=0. 052917720859 nm), energies 
(2i?oo = 219474.63137054 cm'^), and dipole moments (la.u.=2.54174 Debye). 

2. Computational approach 

As in our previous work on alkali dimers [231 [211 [25] or alkali hydrides [26], we 
use an automated procedure to run the CIPSI package (Configuration Interaction by 
Perturbation of a multiconfiguration wave function Selected Iteratively) |27j. The 
approach is based on the gaussian ^-dependent effective core potential (ECP) of Durand 
and Barthelat [281 [29] to represent the ten-electron Mg^"*" ionic core, Gaussian basis sets, 
and ^-dependent effective terms for core polarization (CPP) along the lines of ref. [30] 
revisited by Foucrault et al [31j. The MgH"*" molecule is therefore modeled as an effective 
two-electron system. The molecular orbitals are determined by restricted Hartree-Fock 
single electron calculations, including the CPP [31] providing the potential curves for 
the relevant molecular cations. Then a full configuration interaction (FCI) is performed 
for each relevant molecular symmetry, providing potential curves, and permanent and 
transition dipole moments. 

We will compare our results to those obtained by other various theoretical methods 
falhng in two categories: 
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• Effective core potential methods: Fuentealba and Reyes [32] investigated the ground 
state of alkaline-earth monohydrides ions using a ECP+CPP approach similar to 
ours, and local spin density functional for the valence correlation, instead of a CI. 
In the context of the charge transfer between Mg and H"*" at low-keV energies, 
Dutta et al [17] determined numerous MgH"*" potential curves and couplings, using 
a slightly different ^-dependent gaussian ECP including polarization terms, a basis 
of Slater-type orbitals, and a partial CI. A similar work using a simple exponential 
Hellman-type ECP was previously reported by Numrich and Truhlar [ll]. Finally, 
Garcfa-Madronnal et al [16] performed FCI calculations based on an ECP usually 
referred to as model potential for the Mg"*" ion. 

• All-electron calculations: The most recent all-electron calculations have been 
reported in refs. [33l [T2] , improving the previous similar computations by Monteiro 
et al [19]. The authors employed a multireference configuration interaction 
(MRCI) approach to deliver adiabatic potential curves and dipole moments at 
the multiconfiguration self-consistent field level of the theory (MC-SCF), within 
the commercial MOLPRO package. Prior to this work, Rosmus and Meyer 
[M] studied the MgH"*" ground state with the coupled-electron pair approach 
(CEPA) and CI based on pseudonatural orbitals. Olson and Liu [15j performed 
SCF and CI calculations with Slater-type-function basis set with the ALCHEMY 
package. Spectroscopic constants of the MgH"*" ground state have been calculated 
by Canuto et al [35], using fourth-order-many-body perturbation approach and the 
GAUSSIAN package. 

As a preliminary step, we checked the convergence of our calculations as a function 
of the size of the Gaussian basis sets describing H and Mg"*". 

For hydrogen atom, we considered three basis sets available in the literature. The 
set {H'^"^)} corresponds to the large uncontracted basis {10s5p2(i} that we used for alkali 
hydrides [26j. The basis set {H^^^} ({QsSpld} contracted to {Ispbpld}) has been used 
to study the electronic structure of HC1~ [36], and BaH [37]. The smaller basis set 
{H'^'")} ({7s3p2(i} contracted to {6s3p2(i}) has been used in ref.[38] to investigate the 
structure of alkali hydrides. The deviations for the calculated hydrogen ground state 
binding energy, compared to the experimental one (taken at the non-relativistic value of 
-0.5 a.u.), are Ai, (cm-i)=3.1, 16, 64, for the three basis sets {H^^)}, {H^-^)}, {H^^)}, 
respectively. 

The basis sets for Mg"*" has been chosen (Table [1]) in combination with the 
adjustment of the cut-off radii ps, Pp, Pd, Pf of the CPP term, to match the experimental 
energy of the lowest levels of the ion, or their spin-orbit averaged energy, when 
appropriate (Table [2]). We first used the basis set {Mg^^^^}, proposed in ref. [39] 
{{llsGpld} contracted in {isSpld}) in the study of the structure and the stability 
of Mg+ and Mg++ clusters. The lowest d level of Mg+ cannot be reproduced accurately 
with this basis, which is due to the lack of d orbitals. Therefore we extended the basis 
of ref. [m {{4s2pld} contracted to {3s2pld}) towards the sets {Mg(^)} {{6s4p3d} 
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contracted to {5sAp2d}) and {Mg^*^))} {{7s5p3d2f} contracted to {6s5p2d2f}), which 
are able to reproduce Mg+ experimental binding energies 3s, 3p and 3d exactly, and 
those of the higher excited levels 4s, 4p, 5s at about 1.5% or better. This agreement is 
similar to the one obtained for alkali atoms in our previous works. 

The comparison with the results of refs. flGl [331 [12] can be carried out from the 
point of view of binding energies, or energy differences (i.e. excitation energies, when 
selection rules apply). We first extracted from ref.|T6] atomic binding energies from 
the value of potential energies at R = 30ao. We see that due to their simple form 
of their model potential involving only two parameters, the authors cannot match the 
binding energy of the Mg+ ground level. Moreover they did not include d orbitals in 
their basis, so that their 3d level is indeed ill-defined, as it is visible in Table [2l Their 
lowest energy differences are in better agreement with experiment (again, except the 
one for the 3s — 3d difference) than binding energies, which suggests that their model 
potential leads to a global upward energy shift of the excitation spectrum. Refs. [33l [12] 
report all-electron calculations which cannot provide a precise binding energy for the 
ground state, so that the authors simply shifted the whole atomic spectrum to match 
the experimental ground level energy. Thus they obtained an excitation energy for the 
3s — 3p transition in reasonable agreement (at a 4% level) with experiment. 

Prior to molecular calculations, we first performed a FCI for the neutral Mg atom, 
yielding ground state energies reported in Table [2 We obtain a reasonable agreement 
with the experimental value, comparable to the one obtained for alkali negative ions in 
our previous work [21]. The largest basis set gives a ground state energy lower than 
the experimental one, which may be partly due to negative energy differences obtained 
for the highly-excited levels in Mg+. Our value is in better agreement than the one 
obtained in the other available works. 

3. Potential curve calculations 

The potential curves for the lowest singlet and triplet states of MgH+, computed with 
the combination of the {Mg*^*^))} and {H'^^^)} basis sets, are displayed in Figure [H 
compared those of ref.|16] which are, to our knowledge, the only ones whose numerical 
data is available in the literature. As discussed in the previous section, this comparison 
is presented from the point of view of binding energies. The discrepancies reported in 
Table [2] on atomic energies are easily identified at large distances. It is not surprising 
that well depths are smaller in ref.[T6] than in the present work, due to the limited 
basis set size. This is particularly true for the B^U and the a states, whose potential 
wells found in ref . [16] are hardly visible at the scale of the figure. The splitting of the 
large avoided crossing between the C and the A ^E"*" states, which is responsible for 
charge exchange phenomena at high collision energies [TH [17] , is reduced according to 
our calculations, which could then modify the conclusions for the collision experiments. 

In Figl2l we present the comparison of the four lowest singlet potential curves with 
the recent results of ref. [33l [T2] . In this latter reference, the authors have shifted their 
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Table 1. Exponents (a) and contraction coefficients (b) of the gaussian orbitals of 
the various basis sets presented in the text for Mg+, static polarizability a (in a.u.) of 

Mg"*" and adjusted cut-off radii ps, Pp, pd, Pf (in a.u.) present in the i'-dcpcndcnt CPP 
term. The sets of orbitals which are contracted together are separated by a blank line. 
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Table 2. Energy differences Ecaic — E^xp (in cm^^) between calculated and 
experimental values of the lowest Mg++H(n — 1) asymptotes, and of the Mg(3s^)+H+ 
asymptote, obtained with the various basis sets presented in the text. Experimental 
values are given for binding energies (columns (I) and (II)) and for excitation energies 
(column (III)). Columns (a), (b), (c), hold for the results using basis sets {Mg'"^^)}, 
{Mg(-^)}, and {Mg'*-^)}, respectively. In column (d), energy differences are deduced 
from the values of the potential curves computed in ref.[16] at i? = 30a. u., and from 
the corresponding excitation energies (column (e)). In column (f), energy differences 
are extracted from the excitation energies reported in ref . |33|, [T2] . 

potential curves one by one to the experimental value of their corresponding dissociation 
limit, hoping to minimize the effect of the shift on their computed atomic energies in their 
modeling of the MgH+ photodissociation. Therefore we have drawn all potential curves 
relative to the same origin of energy, in order to emphasize on the differences in the shape 
and depth of the potential curves. The overall agreement between our ECP+CPP-I-FCI 
approach and the all-electron MOLPRO calculations is quite satisfactory, especially for 
the C state, while the well depth of the X, A, and B states is slightly less pronounced. 
The experimental potential curve extracted in ref. [21] from spectroscopic measurements 
concerns only the low part of the well of the X and A states, and we have drawn them 
in order to make them match the bottom of the present curves. The agreement in the 
shape of our potential curves with those of the RKR curves is excellent. 

While complete sets of data for potential curves are quite scarce, several authors 
reported the main spectroscopic constants (equilibrium distance _Re, potential well depth 
Dei ci'iid harmonic constant uj^ for the lowest singlet and triplet molecular states, which 
are summarized in Tables [3] and HI respectively. The main influence of the choice of 
the basis set combination is visible of the potential well depth, and is only weak on the 
harmonic constant, while the equilibrium distance is almost stable. Our recommended 
values would be those yielded by the combination of the {Mg^'-^^)} and {H*^"^))} basis 
sets, while the {H*^^))} set still looks reasonable to use. Our results for the ground 
state are in good agreement with the reported experimental data, except with the one 
of ref. [20]. The situation is less clear for the A state, as we refer to experimental 
data which are quoted in ref. [H], without the original source. Moreover, Balfour [21] 
only provided rough estimates of the depth of the potential curves, as his spectroscopic 
study did not cover the entire potential well. A broad dispersion can be seen in the 
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other theoretical results (apart those of ref s. [TBI [12] already discussed above), and our 
results are in agreement with those of refs. [35| [T9] . Finally, as it could be expected from 
the discussion above, the potential well depth of our triplet states is generally larger 
than the ones available elsewhere [161 [IS]- Note that the c state probably has at most 
a single vibrational level. 
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Table 3. Equilibrium distance Rf, (in a.u.), potential well depth D^, (in cm^^), and 
harmonic constants uj^ for the X^Y^ , ^^1]+, C^i;+, and B ^11 states of MgH+, 
computed in the present work and compared to values available in the literature. 
Calculations with various combinations of basis sets are presented in the uppermost 
rows: (a) {Mg('4))}+{H(^))}, (b) {Mg(^))}+{H(^))}, (c) {Mg(c^))}+{H(^))}, (d) 
{Mg(^))}+{H(-^))}, (e) {Mg(C'))}+{H('^')}. Values displayed in the middle rows are 
from theory, and from experiment in the lowest rows. 

4. Other properties: Transition and permanent dipole moments, 
polarizability 

We have computed the electronic transition and permanent dipole moments, and static 
dipole polarizability, as functions of the internuclear distance, using the various basis 
set combinations defined in the previous section. As in our previous works, we checked 
that the difference among the calculations do not exceed about 1% among each other. 
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Table 4. Equilibrium distance Re (in a.u.), potential well depth (in cm^^), and 
harmonic constants uie for the a'^E^, c^S^, and 6^11 states of MgH+, computed in the 
present work, and compared to values available in the literature. The various displayed 
cases are labeled in the same way than in FiglHl 



in the vicinity of their maximal value. In the figures, we only reported values coming 
from the combination {Mg('^))}+{H(^))}. 

Figure [3] show the dipole moment functions for all the possible transitions between 
the five lowest singlet molecular states. Our results are found in good agreement the 
data from refs. [33t [T2] for the X — A, A — B, and A — C transitions, which is a strong 
argument in favor of the quality of both calculations, which proceed along different lines 
0. Apart from these transitions, the other functions are given for the first time, to our 
knowledge. The A — C transition plays the main role in the photodissociation model of 
ref. [12] . The change of sign in the A — C function around 6 a.u. is due to the change 
of character of the C state due to the avoided crossing with the A state. This function 
rapidly drops down to zero at large distances, as expected for a transition between states 
of the Mg^+H pair and of the Mg+H"*" pair. The maximum of this function around 
5 a.u. coincides with the avoided crossing between the C and D excited states, just like 
in the X — D function. We also note the local change of the A — D function around 
12 a.u. due to the avoided crossing of the D state with an upper excited state. It is 
likely to expect that the efficient excitation of the A state into the D state, even for high- 
lying vibrational levels (i.e. with a wave function extending towards large distances) 
would lead to another photodissociation channel with charge exchange into the C state, 
through predissociation across the avoided crossing between the D and C states. 

The permanent dipole moment of heteronuclear molecular ions is most often tedious 
to derive from standard electronic structure calculations, as its calculated value depends 
on the choice of the origin of coordinates. A straightforward way to overcome this 
difficulty is to compute the potential energy of the molecular ion placed in a perturbative 
static electric field with a constant amplitude and direction. This corresponds to the 
well-known finite- field method suggested by Cohen many years ago |42]. By varying 

X In ref. [12], we changed the sign of the A ~ C function beyond 6 a.u., as the oscillation shown there 
actually corresponds to such a sign change. 
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the strength of the electric field at each distance, the potential energy evolves as a 
parabola, yielding the permanent dipole moment and the static dipole polarizability as 
the linear and quadratic coefficients [25]. If the z axis is chosen along the internuclear 
axis in the molecule-fixed reference frame {x,y,z), they are two independent components 
of the molecular polarizability tensor, i.e., the parallel component a\\ = azz and the 
perpendicular one a± = a^x = ctyy. Two related quantities are usually defined: the 
average polarizability a = {azz+'^o.^x) /"i and the polarizability anisotropy 7 = cy.zz—Oixx- 
Our results for the ground state are shown in Figure HI compared to the very 
few other available theoretical data. The permanent dipole moment is in satisfactory 
agreement with the one of ref.pL2| computed at their highest level of theory, both for 
the position and amplitude of its maximum and for the position of the change of sign. 
However we note that the asymptotic value obtained in ref.jT2] seems to diverge from 
a vanishing value. Both components of the static dipole polarizabilities exhibit an R- 
variation similar to the one obtained for alkali dimers [25], alkali hydrides [26], or for the 
hydrogen molecule [13]. In all these systems, the parallel polarizabilities has a maximum 
at a distance around 1.3 to 1.5 times the equilibrium distance i?e of the system, while the 
perpendicular components always have a smaller magnitude than ay, and monotonically 
increase towards the asymptotic limit. For all the quantities of Figure HJ our results are 
also in agreement with the data reported at the equilibrium distance by Sadlej and Urban 
[B] using polarized basis sets and many-body perturbation and coupled-cluster theories 
combined with finite field method. At large distances, the polarizability components 
converge toward the sum of the atomic values of the two constituents. We found 
an asymptotic value (39.9 a.u.) which is slightly smaller than the sum of the atomic 
polarizabilities for H (4.5 a.u. ^) and for Mg+ (38.84 a.u.) [44J. 

5. Conclusion 

In the present paper, we performed an accurate investigation of the electronic structure 
of the MgH"*" ion, including potential curves for singlet and triplet states, transition and 
permanent dipole moments, and static polarizabilities, as functions of the internuclear 
distances. Our work extends the recent similar study of ref.[T2] using an all-electron 
approach based on the MOLPRO package, restricted to the lowest singlet states. While 
minor differences occurred, both approaches are in good agreement, assessing their 
validity and their good accuracy. By computing higher electronic states, our study 
suggests that another MgH"*" photodissociation path could be worthwhile to try in the 
framework of the trapped Mg+ experiments involving collisions with H2, resulting into 
the formation of trapped MgH"*" molecular ions. Our study represents the necessary 
initial step towards the description of the quantum dynamics of the Mg"'"(3p)+H2 
reactivity, which has been recently probed to exhibit a strong isotopic effect [13]. 
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R (a.u.) R (a.u.) 



Figure 1. (a) Singlet and (b) triplet potential curves MgH+ computed in the 
present work (full lines: S"*" symmetry; dotted-dashed lines: 11 symmetry) compared 
to those of ref.[T6] (dashed lines: E+ symmetry; dotted lines: 11 symmetry). The 
lowest states are labeled by their standard spectroscopic notation: X, A, B,C, D for 
(l)iS+(3s), (2)il]+(3p), (l)in(3p), (3)iS+(3s2), (1) (4s), respectively, and a, 6, c 
for (l)3s+(3s), (l)3n(3p), (2)3l]+(3p), respectively 
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R(a.u.) 

Figure 2. Potential curves of the four lowest singlet states of MgH+, using the same 
spectroscopic notation than in FiglT] Full lines: present calculations. Dashed lines: 
ref.[33l [12], at the highest level of theory reported there. All curves are shifted to 
the same origin, for better visibility. The experimental RKR curves for the X and A 
states from ref. 2l] are drawn with circles. They are shifted to the bottom of the wells 
calculated in the present work. 




Figure 3. Transition dipole moments for selected states of MgH+. Transitions are 
labeled using the same spectroscopic notation than in FiglTJ Lines without symbol 
correspond to the present calculations, while those marked with circles are from 
ref. [33l [12] , computed at their highest level of theory . 




Figure 4. (a) Permanent dipole moment of the MgH+ ground state computed in tlie 
present woric (full line) compared to the one of ref.ref.[331 [H] (dashed line)computed 
at their highest level of theory, (b) Parallel (full line) and perpendicular (dashed 
line) dipole polarizability of the MgH+ ground state computed in the present work. 
Theoretical values reported in ref.[44 at the equilibrium distance (3.12ao) are indicated 
by stars and triangle. 



